arXiv:1507.01607v2 [hep-ph] 26 Sep 2015 
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Recent experimental data on the T(45') —>■ T{lS)ri and T(4S) —>■ hb(lP)r] processes seem 
to contradict the naive expectation that hadronic transitions with spin-flipping terms should be 
suppressed with respect those without spin-flip. We analyze these transitions using the QCD 
Multipole Expansion (QCDME) approach and within a constituent quark model framework that has 
been applied successfully to the heavy-quark sectors during the last years. The QCDME formalism 
requires the computation of hybrid intermediate states which has been performed in a natural, 
parameter-free extension of our constituent quark model based on the Quark Confining String 
(QCS) scheme. We show that i) the Ml-Ml contribution in the decay rate of the T(4S) ^ T(1S')?7 is 
important and its suppression until now is not justified; ii) the role played by the L = 0 hybrid states, 
which enter in the calculation of the Ml-Ml contribution, explains the observed enhancement in the 
T(4S) —>■ T(15')77 decay width; and hi) the anomalously large decay rate of the T(45') —>■ hb{lP)ri 
transition has the same physical origin. 
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I. INTRODUCTION 


The general way of referring to an hadronic transition 
is [1] 


<Pi^^F+h, (1) 

where <i>/ and <1)^^ stand, respectively, for the initial 
and final states of heavy quarkonium. The emitted 
light hadron(s), h, are kinematically dominated by single 
particle (tt^, r],uj, ...) or two particle (27r, 2K, ...) states. 

Hadronic transitions are important decay modes for 
low-lying heavy quarkonium states. For instance, the 
first observed hadronic transition'0(25') —>■ J/^tttt [2] has 
a branching fraction reported by the Particle Data Group 
(PDG) of (52.58 ± 0.43)% [3]. Moreover, during the 
last years, hadronic transitions between heavy quarkonia 
have led to a remarkable series of discoveries helping 
either to establish new conventional heavy quarkonium 
states, like the hb{lP) and hb{2P) observed in the two- 
pion decay of the T(55) [4] state, or to extract relevant 
information of the so-called “XYZ” states, like in the 
cases of X(3872) [5-7] and X(4260) [8]. 

The BaBar Collaboration has presented a systematic 
study of hadronic transitions between T{mS) (m = 
4,3,2) and T{nS) (n = 2,1) states, reporting the first 
observation of the T(45) —>■ T{lS)r] [9]. The measured 
branching fraction for this decay is 

H(T(45) -> T(15)77) = (1.96 ± 0.06 ± 0.09) x 10-^ (2) 
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which is puzzling larger than the branching fraction for 
T(45) T(15)7r+7r-, (0.800 ± 0.064 ± 0.027) x 10-^ 

with a ratio between them of 


R,[T(45)] 


r(T(45) ^ T(15)r;) 
r(T(45) -> T(15)7r+7r-) 


2.41±0.40±0.12. 


. . 

The comparison of the T(mS) —>■ T(nS)r] transitions, 
which are spin-flipping, and the corresponding tt+tt” 
ones, that do no require the 5-quark spin to flip, is 
particularly interesting because we naively expect a 
suppression of the transitions with spin-flipping terms. 
This is observed for the two vector bottomonium states 
that are below open 5-flavored threshold [3] : 


R^[T(25)] = (1.63 ± 0.23) x 10■^ 
i?^[T(35)] < 2.29 X 10-3, 


but also in the charmonium sector with i?^[0(25)] = 
(9.75 ± 0.17) X 10“^ [3]. Note that even for the state 
which is just close above the c-flavored threshold, and 
thus it should present a similar role in the charmonium 
sector than the T(45) in the bottomonium one, the PDG 
reports a value of i?,,[0(3770)] = 0.47 ± 0.22. 

Further insight into the anomalously large T(45) —>■ 
T(15)77 decay rate can be gained by searching for the 
transition T(45) —>■ hb{lP)r/ because it is, in principle, 
dominated by similar spin-flipping contributions. The 
Belle Collaboration has very recently measured for the 
first time the branching fraction H(T(45) —>■ hb(\P)r]) = 
(2.18 ±0.11 ±0.18) X 10-3 [10]. This branching fraction 
provides a partial decay width of (44.69±6.95) keV when 
combined with the total decay rate reported by PDG [3]. 
This value is again unexpectedly large. 

The anomalous hadronic decay widths can be due to 
several mechanisms: Contribution of hadron loops [11- 
13]; four-quark components in the quarkonium wave 
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functions [14]; internal loop radiation [15, 16]; or, as we 
have pointed out in Ref. [17] and we will see herein, the 
existence of hybrid mesons with a mass near the one of 
the decaying resonance^. 

The standard theoretical approach to study hadronic 
transitions is QCD Multipole Expansion (QCDME) [19- 
24] . Tung-Mow Yan was the first one to present a gauge- 
invariant formulation within this framework [25]. Many 
details about QCDME in the context of the Kuang-Yan 
model can be found, for instance, in Refs. [26, 27]. The 
interested reader is also referred to the recent review [28] 
of Yu-Ping Kuang. 

This approach describes the hadronic transition as 
a two-step process in which the heavy quark system 
initially emits, at least, two gluons that subsequently 
hadronize into light hadrons. After the emission of the 
first gluon and before the emission of the second one, 
there exists a propagating intermediate state where the 
QQ pair together with the gluon forms a hybrid state. 

Other possibility is to work in a local approxima¬ 
tion [24] which does not require the description of hybrid 
states but it is strictly valid only in the limit of infinite 
heavy quark mass. 

The width of a hadronic transition in QCDME depends 
critically on the position in the spectrum of the hybrid 
states, therefore it is important to describe consistently 
the heavy quarkonia and the hybrids using as few 
parameters as possible. A description of hybrid mesons 
is difficult from first principles of QCD [29, 30] and 
one is generally forced to use models: the flux-tube 
model [31, 32], constituent gluons [33], Coulomb gauge 
QCD [34], quark confining string model (QCS) [35-38] 
or QCD string model [39]. 

The hadronization vertex is independent of the prop¬ 
erties of heavy quarkonium and hybrid states. Usually, 
this kind of vertices are calculated using low energy the¬ 
orems [40] . This procedure limits the predictive power of 
the approach since for each new configuration of initial 
and final bound states one needs new adjustable param¬ 
eters. Moreover, for the interesting case of the hadronic 
transitions involving the 77 meson one must discard cer¬ 
tain operators involving the gluon fields. This last point 
is relevant for the discussion of this work. 

One important feature of the QCDME approach is 
that the hadronic transitions are classified in a series 
of chromoelectric and chromomagnetic multipoles. The 
decay rate of the T(45') —> T( 1 S ')?7 process has two 
leading multipole gluon emission terms: Ml-Ml and El- 
M2, whereas the leading order term in the T(45') —>■ 
T(lS') 7 r 7 r transition is an El-El contribution. The Ml- 
Ml term has been usually neglected. However, the 
hybrid mesons involved in the calculation of the Ml-Ml 
amplitude are different than those involved in the El- 
M2/E1-E1 contributions. Therefore, we consider worthy 


to estimate in the T(4S') —)■ T( 1 S ')77 process the Ml- 
Ml term and thus the effect produced in the amplitude 
by the hybrid meson spectrum before resorting to more 
sophisticated or exotic mechanisms. 

To do this, we will use the QCDME within the 
same theoretical formalism presented in Ref. [17]. This 
formalism has explained successfully some puzzles of 
the two-pion hadronic transitions in the charmonium 
sector. The branching ratios i?^[T(n5')] with n = 2, 3, 4 
will be presented herein^. An important feature is 
that our model for hybrid mesons [17] is a natural, 
parameter-free extension of a constituent quark model 
(CQM) [41] (for reviews on the CQM, see Refs. [42, 43]) 
that describes quite well hadron phenomenology and 
hadronic reactions [44-46]. Eurthermore, the CQM has 
been recently applied to mesons containing heavy quarks 
with a remarkable success, describing a wide range of 
physical observables which concern spectrum [47-49], 
strong reactions [50-52] and weak decays [53-55]. 

This manuscript is arranged as follows. In Sec. 2 
we introduce the constituent quark model pointing out 
only those features which are relevant for this work. 
Section 3 is dedicated to explain the parameter-free 
extension of our quark model to describe hybrid mesons. 
Section 4 shows the QCDME formulation of the hadronic 
transitions highlighting its most relevant features. Our 
results are provided in Sec. 5. We finish giving some 
conclusions in Sec. 6 . 


II. CONSTITUENT QUARK MODEL 

The description of the meson spectra is based on the 
constituent quark model (CQM) proposed by Vijande 
et al. in Ref. [41]. One must solve the Schrbdinger 
equation with a quark-antiquark potential whose main 
pieces are: i) the Goldstone-boson exchanges between 
dressed constituent quarks which is a consequence of 
the dynamical chiral symmetry breaking of QCD, ii) the 
perturbative one-gluon fluctuations around the instanton 
vacuum, and iii) a phenomenological confining potential 
which reflects the empirical fact that quarks and gluons 
have never seen as isolated particles. Note that in the 
heavy quark sector chiral symmetry is explicitly broken 
and thus Goldstone-boson exchanges do not appear. 

Further details about the CQM and the fine-tuned 
model parameters can be found in Refs. [41, 47, 56]. 
Here we want to explain in more detail our confinement 
potential because its screened linear shape is a particular 
feature of the model. It is well known that multigluon 
exchanges produce an attractive linearly rising potential 
proportional to the distance between infinite heavy 
quarks. However, sea quarks are also important 


^ A discussion about the goodness of the approach applied to the 
calculation of the hadronic decays treated herein can be found 
in Appendix A. 


^ A similar observation has been also done in Ref. [18]. 
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State 

Mass (CQM) Mass (Cornell) 
(MeV) (MeV) 

Mass (Exp.) 
(MeV) 

T(1S) 

9502 

9460 

9460.30 ± 0.26 

T(2S') 

10015 

10050 

10023.26 ±0.31 

r{3S) 

10349 

10400 

10355.2 ±0.5 

T(4S) 

10607 

10670 

10579.4 ± 1.2 


Table I. Masses, in MeV, of the S-wsve vector bottomonium 
states up to n = 4 predicted by our CQM and by the Cornell 
model [58, 59]. Experimental masses are taken from PDG [3]. 
The mass of the T(1S) has been fitted in the Cornell potential. 


ingredients of the strong interaction dynamics that 
contribute to the screening of the rising potential at low 
momenta and eventually to the breaking of the quark- 
antiquark binding string [57]. Our model try to mimic 
this behavior using the following expression 

VcoMr) = [-ac(l - + A] (A^ • A|), (5) 

where Uc and are parameters, r is the interquark 
distance and A°^-j are SU{3) color matrices. At short 
distances this potential presents a linear behavior with 
an effective confinement strength a = —acfJ.ci^q ’ 
while it becomes constant at large distances showing a 
threshold defined by fthr = [—Oc + ’ ^q)- No qq 

bound states can be found for energies higher than this 
threshold. 

The screened linear potential is a key feature to 
reproduce the degeneracy pattern observed for the higher 
excited states of light mesons [56]. As we assume that 
confining interaction is flavor independent, this affects 
also to the different quark sectors and, in particular, 
the bottomonium spectrum. Table I shows the S- 
wave vector bottomonium states up to n = 4 predicted 
by our quark model. We compare masses with the 
standard Cornell model [58, 59] and the experimental 
data reported by PDG [3]. Both models reproduce quite 
well the experimental data. In the Cornell model the 
ground state is fitted to the experimental figure while in 
our case the model parameters are fitted to all meson 
sectors. One can see that our quark model is able to 
reproduce in better agreement the excited states which 
is a particular feature of the screening of the linear 
confinement potential. 


III. A MODEL FOR HYBRIDS 

From the generic properties of QCD, we might expect 
to have states in which the gluonic field itself is excited 
and carries quantum numbers. A bound-state is 

called glueball when any valence quark content is absent, 
the addition of a constituent quark-antiquark pair to an 
excited gluonic field gives rise to what is called an hybrid 
meson. The gluonic quantum numbers couple to those 
of the qq pair. This coupling may give rise to so-called 


exotic mesons, but also can produce hybrid mesons 
with natural quantum numbers. We are interested on 
describing the last ones because they are involved in the 
calculation of hadronic transitions within the QCDME 
approach. 

Ab-initio QCD calculations of the hybrid (even con¬ 
ventional) bottomonium states are particularly difficult 
because the large mass of the 6-quark. Therefore, the 
only way up to now to describe hybrid mesons in the 
bottomonium sector is through models. An extension 
of the quark model described above to include hybrid 
states has been presented in Ref. [17]. This extension 
is inspired on the Buchmuller-Tye quark-confining string 
(QCS) model [35-38] which assumes that the meson is 
composed of a quark and antiquark linked by an appro¬ 
priate color electric flux line: the string. Gluon excitation 
effects are described by the vibration of the string. These 
vibrational modes provide new states beyond the naive 
meson picture and are interpreted as hybrid mesons. 

The coupled equations that describe the dynamics 
of the string, quark and antiquark sectors are highly 
nonlinear so that there is no hope of solving them 
completely. Then, to introduce the vibrational modes, 
we use the following approximation scheme. First, we 
solve the string Hamiltonian via the Bohr-Oppenheimer 
method to obtain the vibrational energies as a function 
of the interquark distance [36] 


Vn{r) = cr(r)r 



2mT 

a{r) [(r — 2dY + dd^] 


1/2 


( 6 ) 


Note that n = 0 gives Vb(r) = (j{r)r where cr(r) = 
(16/3) Oc [(1 — e“'^‘=’')/r] attending to Eq. (5). The 
parameter d is the correction due to the finite heavy 
quark mass 


d{mQ,r, a,n) 


ar^an 


4(2mQ -I- crra„) ’ 


(7) 


where is related with the shape of the vibrating string 
and can take the values 1 < q;„ < \/2. 

Second, the vibrational potential is inserted into the 
meson equation as an effective potential 


yhyh{r) = VbGE(?') + VcoN{r) + [Vn{r) - a{r)r] , (8) 

where VbGE('c) + Vcon('c) is the naive quark-antiquark 
potential in the heavy quark sector with VbGE(?') = 
— 4as/3r and Vcon(?’) given in Eq. (5). Vn(r) is the 
vibrational potential calculated above. We must subtract 
the term a{r)r because it appears twice, one in Vcon(?') 
and the other one in Vn{r). 

We have arrived to a description of the hybrid mesons 
in the heavy quark sector that does not include new 
parameters besides those of the original quark modeP. 


® The variation of the parameter within its range [1, \/2] 
modifies around 30 MeV the mass of a hybrid state in the 
bottomonium sector, we have chosen the mean value a„ = Vl-S. 
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K / L 

0 

1 

1 

10571 

10785 

2 

10857 

10999 

3 

11063 

11175 

4 

11232 

11325 

5 

11374 

11452 

6 

11496 

11562 

7 

11600 

11657 

8 

11690 

11738 

9 

11766 

11807 

10 

11831 

11866 

11 

11885 

11913 

12 

11927 

- 


Khr = 11943 


Table II. Masses, in MeV, of those hybrid states in the bot- 
tomonium sector that participate in the hadronic transitions 
calculated in this paper. K and L are the hybrid meson quan¬ 
tum numbers. 





Figure 1. A hadronic transition as a two-step process: (1) 
emission of gluons from heavy quarks (MGE), and (2) the 
conversion of gluons into light hadrons (H). 


In that sense, our calculation of the hybrid states is 
parameter-free. Another important feature of our hybrid 
model is that, just like the naive quark model, the hybrid 
potential has a threshold from which no more bound 
states can be found and so we have a finite number of 
hybrid states in the spectrum. 

Table II shows our theoretical prediction for those 
hybrid states in the bottomonium sector that will 
participate in the hadronic transitions in which we are 
interested. They are classified by the angular momentum 
L and radial excitation K. It is important to realize here 
that we predict a hybrid meson with quantum numbers 
\KL) = |10) which is very close in mass, 10571 MeV, to 
the T(4,S) state, 10607 MeV. As we will see later, this 
feature have important consequences in the calculation 
of the T(4S') hadronic decay rates. 


IV. HADRONIC DECAY RATES 


The Hamiltonian for a heavy QQ system in QCDME 
is given by [25]: 

zyeff _ 4- _L iQl 

^QCD — -"QCD 4" ^QCD '"QCD: 1^1 

where represents the sum of the kinetic and 

potential energies of the heavy quarks, is related 

with the color charge of the QQ system (which is zero for 
color singlets), and 77 qcd couples color singlets to octet 
QQ states. Therefore, the hadronic transitions between 
eigenstates |4>/) and of 77 qcd are at least second 


(‘ 2 '] 

order in Tfgcjj and the leading term is given by 




_ 

P 4- if) — QCD 

(2) |ivr\ _ 


= 


^ (4>ph| \KL) {KL \ilgo 1$/), 

( 10 ) 


KL 


where \KL) with associated energies EpL are intermedi¬ 
ate states after the emission of the first gluon and before 
the emission of the second one. They are the hybrid 
mesons described in the former Section. 

A connection is made to the physical process in Eq. (1) 
by assuming that the hadronic transition amplitude 
always splits into two factors (see Fig. 1). The first one 
concerns the multipole gluon emission (MGE) from the 
heavy quarks and the second one is an hadronization (H) 
process describing the conversion of the emitted gluons 
into light hadron(s). 

The MGE vertex involves the wave functions and 
energies of the initial and hnal quarkonium states as well 
as those of the intermediate hybrid mesons. All these 
quantities are calculated within our model and enter in 
the hadronic transition rate with integrals of the type^: 


rLPjPp 

IIP 


=E 


1 


K 


Mi - Mkl 


J dr Rp{r)RKL (r) 


dr' 




RKL{'r')Ri{r') 


( 11 ) 


^ A detailed description of the computation of the decay rates 
in the single-channel approach of QCDME for the hadronic 
transitions can be found in Ref. [26]. See also Ref. [28] for an 
updated review and Ref. [17] for a calculation of two-pion spin- 
nonflip hadronic transitions within our model. 
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where Ri{r) and Rpir) are, respectively, the radial wave 
functions of the initial and final states. RKbij) is 
the radial wave function of the intermediate vibrational 
states \KL). The mass of the decaying meson is M/, 
whereas the ones corresponding to the hybrid states are 
Mkl- 

The H vertex is at the scale of the light hadron(s) and 
is independent of the properties of the heavy quarkonia 
and hybrid states. There are two ways of calculating 
the matrix elements associated with the H factor: (HI) 
using PCAC and soft pion techniques [25, 40] or (H2) 
approximating the hadronic transition rates by 2 -gluon 
emission rates [26], for instance, r(<l>/ —^ ^F'i) — 
r(<l >7 —>• ^F{gg)o-) where the two gluons are projected 
into = 0 “ to simulate the g meson. 

Certainly, the H2-approach does not take in very 
detail the conversion of gluonic field(s) into a single g 
meson. This can be view as a crude approximation and 
one should expect large uncertainties in its predictions. 
However, it is the only formalism that allows us to treat 
consistently the Ml-Ml amplitude. One could estimate 
the uncertainties introduced by the H2-approach using 
the g fragmentation function for the gluons at the scale 
of interest. This universal nonperturbative object has 
been calculated, e.g., in Ref. [60] at NLO accuracy and 
at the scale g = 1 GeV. The range of values significant 
for the processes studied here is between 0.25 and 0.50. 
Therefore, the results predicted by the H2-approach 
could have a systematic uncertainty of about 50% but 
not orders of magnitude that can change drastically our 
conclusions. 

The HI and H2 approaches involve unknown coeffi¬ 
cients. However, the difference between the two is that 
the Hl-approach needs, at least, one coefficient for each 
multipole matrix element (we will denote these coeffi¬ 
cients by Ci’s), whereas in the H2-approach all these 
matrix elements are written in function of only two pa¬ 
rameters {gs and gm)- 

We have mentioned in the Introduction that the 
leading multipoles of an g transition between spin-triplet 
S'-wave states are Ml-Ml and E1-M2. Therefore, the 
matrix element is given schematically by 

A4(^«S'i —>■ ^51 + g) = M-mimi + M.eim 2 - ( 12 ) 

The most extended version of the QCDME formalism, 
the Hl-approach, only considers the E1-M2 multipole 
gluon emission because it is difficult to find a relation 
to fix the Ml-Ml corresponding parameter, whereas the 
H2-approach considers both of them. There is no reason 
to neglect the Ml-Ml contribution and this will be a 
key feature in order to explain the large value of the 
T(4S') —>• T(lS')r 7 decay rate. 

For completeness and because it is important for 
the discussion in the next Section, we show here the 
expressions in both approaches for the —>■ -|- tttt 


transition: 


^ TTTT 

pi/2 
^ TTTT 


(^\ (Alj — MfY I j111|2 
V 2 y 18907r3 ' ’ 


(13) 

(14) 


and for the 35*;^ —>• 35 ^ _|_ ^ transition: 




27 Miml 


I IF 


t^H2 _ 

1 T] 


127r3 


9m \ 

f gEgM \ ‘ 

\ Smg / 127r3 


1 [Ml MfY I J.00Q12 


140 


I/: 


IF 


6804 


(15) 


1 {Mj-Mf? 

\JlF \ I 


(16) 


where the factor G is a phase-space integral defined in 
Eq. (2.4) of Ref. [26], g is the momentum of g and 
has been defined in Eq. (11). 


V. RESULTS 


The transitions T(2S') —>■ T(15')7r+7r“ and T(25') —>■ 
T( 1 S ')?7 help us to fix our unknown coefficients. The 
experimental data [3] 

r(T(25') ^ T(15')7r+7r-) = 5.71 ± 0.48keV, (17) 


fixes the constant Gi in the Hl-approach. Once we get 
Cl the coefficient gs is determined through Eqs. (13) 
and (14). Our values of these two constants are^: 


Gi = 7.69 X 10■^ 
gE = 2.17, 


(18) 


which compare well, for instance, with the ones used in 
the Kuang-Yan model [26]. 

Now, the constants G 3 and gM are fixed to the 
experimental figure of the T(25') —>■ T(1S')?7 decay 
rate [3] 

r(T(2S’) ^ T(1S')77) = (9.27 ± 1.49) x lO'^ keV, (19) 


and we obtain: 


C 3 = 2.96 X 10+® 
gM = 5-70 


( 20 ) 


which are also compatible with the values reported in 
Ref. [26] using the Kuang-Yan model. 


3 Note that the value reported here for Ci differs to the one used 
in our previous work [17]: 9.69 X we decided to fit here 

the bottomonium case instead of the charmonium one in order 
to eliminate sources of uncertainty. One can realize that the 
difference is small (~ 10%) and it is not going to change our 
conclusions. 
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Process Hl-approach H2-approach Experiment [3] 

T(25) ^ T(15)7r+7r- (keV) 5?rL 571 5.71 ±0.48 

T(25) ^ T(15)r; (keV) 9.27 x 10^^ 9.27 x 10^^ (9.27 ± 1.49) x 10^=* 

ii,[T(25)] 1.62x10^^ 1.62x10^^ (1.64 ± 0.25) x 10^^ 

T(35) ^ T(15)7r+7r- (keV) L18 080 0.89 ± 0.08 

T(35) ^ T(15)r; (keV) 7.59 x 10^^ 20.58 x 10^^ < 2.03 x 10^^ 

ii,[T(35)] 6.43x10^^ 25.7x10^^ <2.29x10^^ 

T(45) ^ T(15)7r+7r- (keV) 4.01 2.54 1.66 ±0.24 

T(45) T(15)r; (keV) 12.60 x 10^^ 6.05 4.02 ± 0.76 

fl,[T(45)]3.14 X 10^=^ 2.38_2.42 ± 0.39 


Table III. The decay widths, in keV, of the T(3S) and T(4S') hadronic transitions into T(lS)7r"'"7r and T(15')r; channels. The 
ratio J?^[T(n5)] = r(T(n5') — 1 T(15')r7)/r(T(n5') —t T(lS)7r’^7r~) is also given. The experimental data is taken from PDG [3]. 


Table III shows our theoretical results for the hadronic 
transitions of the T(3S') and T(4S') into T(lS')7r+7r“ and 
T(15)77 channels. The ratio I?^[T(nS')] = r(T(n5') —>■ 
T(15')?7)/r(T(nS') ^ T(lS')7r+7r-) for n = 3, 4 is also 
given. The most remarkable feature of this table is the 
result we obtain for the i?^[T(4S')]. While the ratio 
predicted within the Hl-approach is of the order of 10“^, 
the calculated value using the H2-approach is of the 
order of unity and agrees nicely with the experimental 
measurement. 

Our result has a natural explanation. The rate of the 
T(45') hadronic decay into the T{lS)ri final state has 
two terms. One is a Ml-Ml contribution which involves 
the hybrid states of T = 0 through the term (see 
Eq. (16)), whereas the other one is a E1-M2 contribution 
which involves the L = 1 hybrid mesons through the term 
(see again Eq. (16)). Therefore, different hybrid 
intermediate states enter in the calculation of the Ml- 
Ml and E1-M2 terms. Our result indicates that the Ml- 
Ml contribution to the hadronic decay rate can be (very) 
important because the role played by the L = 0 hybrid 
spectrum. 

The L = 0 hybrid spectrum should be lower in energy 
than the L = 1 spectrum. This is a general feature, 
but a model-dependent result is given in Table II. The 
masses predicted by our model for the ground states of 
hybrid mesons with L = 0 and L = 1 are 10.6 GeV 
and 10.8 GeV, respectively. The T(4S') is very close in 
mass to the ground state of hybrid mesons with L = 0, 
and thus the mass denominator Mj — Mkl in Eq. (11) 
leads to a large enhancement of the T(4S') —>■ T(15')?7 
decay rate through the Ml-Ml term. The mass splitting 
between the T (dS") and the ground state of hybrid mesons 
with L = 1 is around 0.2 GeV, which is enough to 
predict r(T(45') —>• T(15')77) = 12.92 x 10“^keV when 
we consider only the E1-M2 term in the H2-approach. 
Note that this number is very close to the one obtained 
using the Hl-approach. 

Table HI also shows the T(3S') —>■ T(lS')7r+7r“ and 
T(35) — T{lS)'q decay rates using HI- and H2- 
approaches. The theoretical widths of the T(3S') —>■ 
T(15')7r“''7r“ process are compatible with the experimen¬ 
tal data, being that of the H2-approach in better agree¬ 


ment with experiment. Only an upper limit of the de¬ 
cay rate for the T(35') —>■ T{lS)r] process is reported 
by PDG [3]. Our predicted values are higher than 
this upper limit by a factor ^ 3 for the Hl-approach 
and a factor ^ 10 for the H2-approach. The factor is 
larger in the H2-approach because there are two con¬ 
tributions to the decay rate, Ml-Ml and E1-M2. Tak¬ 
ing into account only the E1-M2 contribution we obtain 
r(T(3S') —> T(15')77) = 2.59 x 10“^ which is compati¬ 
ble with the experimental upper limit but still slightly 
higher. We encourage experimentalists to determine this 
decay rate with better precision in order to clarify the 
situation. 

A couple of comments are necessary here related with 
the existence in the literature of alternative explanations 
to the anomalously large T(45') —>■ T(1S')?7 decay rate. 
The enhancement of the branching ratio i?,,[T(45')] 
could be attributed to neglect the effect of heavy-meson 
loops (see the related discussion for the charmonium 
sector in Refs. [1, 13]). A nonrelativistic effective field 
theory (NREFT) was introduced in Ref. [12] with the 
goal of determining the effect of heavy-meson loops on 
the hadronic transitions between heavy quarkonia with 
controlled uncertainty. The expansion parameter is the 
velocity of the heavy mesons in the intermediate state: 
■cioop- The NREFT determines that a typical hadronic 
transition via heavy-meson loop scales as 

Doop/(Doop)^ X vertex factor, (21) 

and, in the case of i?,,[T(4S')], the vertex factor accounts 
for the transition between two S'-wave bottomonia which 
takes place through a BB loop via a P-wave vertex. 
Therefore, the vertex factor is proportional to 'cfoop and 
the heavy-meson loop contribution scales as order nioop 
with 


'^loop 


y|,„[T( 45 ) - 2 mm ^ , 26 x 10-, 

m[B] 


( 22 ) 


Then, heavy-meson loops cannot explain the order of 
magnitude measured in the branching ratio i?,,[T(4S')]. 

An alternative explanation for the large decay width 
of the T(45') —>■ T(15')?7 transition has been reported 
in Ref. [16]. Based on the Field Correlator Method 
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(FCM) [15], the authors of [16] assume that the rj 
transition between two heavy quarkonia proceeds via 
intermediate states of BB, BB*, etc. with the rj meson 
emitted simultaneously at vertices. The dynamics of 
FCM is different than that of the QCDME and thus it can 
produce a different outcome for the decay widths. The 
results within FCM approach can be found in Table 3 
of Ref. [16]. The order of magnitude is reproduced for 
the r(T(4S') —^ T(lS')r 7 ) but a very large value, in 
strong disagreement with experiment, is also obtained 
for the decay rate of the T(3S') —?► T(1S')?7 transition. In 
order to alleviate this discrepancy, the authors modify 
their model parameters by (10 — 15)%. Surprisingly, 
this keeps unmodified all partial decay widths except 
the r(T(3S') —>■ T(1S')?7) that is reduced by 3 orders 
of magnitude due to cancellations (see again Table 3 of 
Ref. [16]). 

A way to test our picture is computing other hadronic 
transitions of the T(4S') meson in which the L = 0 
hybrids are involved. One example is the T(4iS') —>■ 
hi,{lP)r] transition that has recently been observed by the 
Belle Collaboration [10] and whose measured branching 
fraction appears to be anomalously large. The expression 
for the decay rate of the T{mS) —^ hb{nP)ri hadronic 
transitions can be extracted from Eq. (56) of Ref. [28] 
and reads as 


r(T(mS') ^ hb{nP)rj) = 

_ gh MF\q\ ( 47r 


1144m^ M/ 




I/, 


001 

IF 


110|2 


//F 


(23) 


As one can see in Eq. (23), both L = 0 and L = 1 
hybrid states are involved in this kind of transitions 
through the factors and fjjP, respectively. We 
predict 


F(T(4S')h6(lP)77) = 37.89keV, (24) 


which is in good agreement with the experimental figure, 
(44.69 ± 6.95) keV. The enhancement of this decay width 
has the same physical origin than the one found in the 
decay rate of the T(45') — T{lS)r] process: we are 
predicting a hybrid bottomonium meson which is very 
close in mass to the T(4S'). The location of the hybrid 
mesons in the spectrum is a model-dependent result 
but our prediction seems to explain two uncorrelated 
processes, T(4S') —T( 1 S ')?7 and T(4S') —>■ hb{lP)'rj, 
even a third one considering our results in Ref. [17]. It 
is worth to remark again that the construction of our 
hybrid model does not add any new parameter of a quark 
model which explains a very large number of hadron 
phenomenology. 


VI. CONCLUSIONS 

Recent experiments on the T(4S') —>■ T(1S')?7 and 
T(45') —>■ hb(lP)r] transitions have pointed out their 


anomalously large decay rates. This seems to contradict 
the naive expectation that hadronic transitions with spin¬ 
flipping terms should be suppressed with respect those 
that do not have these terms. 

We have studied these transitions within the theoret¬ 
ical framework of a constituent quark model that has 
been applied successfully to a wide range of hadron ob¬ 
servables. In particular, this model has been used in 
the last few years to study spectra, strong reactions and 
weak decays in the heavy quark sectors. Therefore, we 
consider that the quark model parameters are very well 
constrained for the study of hadronic transitions between 
heavy quarkonia. 

The calculation of the hadronic decay rates has been 
performed using the QCDME approach. This formalism 
requires the computation of a hybrid meson spectrum. 
We have calculated the hybrid states using a natural, 
parameter-free extension of our quark model based on 
the Quark Confining String scheme. 

The rate of the T(4S') —>■ T(15')77 hadronic decay has 
two leading multipole gluon emission terms: Ml-Ml and 
EI-M2. The Ml-Ml term has been usually neglected. 
We have shown that the hybrid mesons involved in the 
calculation of the Ml-Ml and E1-M2 contributions are 
different. The Ml-Ml contribution involves the hybrid 
states of L = 0, whereas the E1-M2 contribution involves 
the L = 1 hybrid mesons. The T(4iS') is very close in 
mass to our predicted ground state of hybrid mesons 
with L = 0, and this leads to a large enhancement of the 
T(4S') —>■ T{lS)ri decay rate through the Ml-Ml term. 

A way to test our prediction against others is looking 
for hadronic transitions that involve the L = 0 hybrids 
and see if similar enhancements have been observed 
experimentally. An example is the anomalously large 
decay rate of the T(4S') —>■ hb{lP)ri process measured 
recently by the Belle Collaboration. We have shown 
here that this enhancement has the same physical origin 
than the one found in the decay T(4S') — T{lS)r]. We 
must admit that the location of the hybrid mesons in the 
spectrum is a mo del-dependent result but our prediction 
seems to explain the two processes above and even a third 
one considering our results in Ref. [17]. 
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Appendix A: GOODNESS OF THE QCDME 
APPROACH 


The QCDME approach is based on performing an 
expansion of the gluon field t) in Taylor series of 

{x — X) at the center of mass position A®: 

Ag(f, t) = Ag(A, t)-{x-X)- E“(A, t) + ..., 

1 ... (Al) 

= — (x-X) X B “( A , t ) + ..., 


where if® and 5“ are color-electric and color-magnetic 
fields, respectively. 

The Hamiltonian formulation is more convenient when 
one wants to follow a nonrelativistic formalism. The cor¬ 
responding Hamiltonian derived from the above formu- 
lation is given in Eq. (9) where, more explicitly, 
is 

Efgo = -da ■ E“(A, t)-ma- B^{x, t) + ..., (A2) 

with 


da = gE J{x - X)^^{x,t)^'^{x,t)d^x, 

ma = ^ J{x- X) x^\x,t)j^'^{x,t)d^x, 


the color-electric dipole moment (El) and the color 
magnetic dipole moment (Ml) of the QQ system, 
respectively. Higher-order terms (not shown) give rise to 
higher-order electric (E2, E3, ...) and magnetic moments 
(M2, ...). 

(2l 

Since couples color-singlet to octet quark- 

antiquark states, the hadronic transitions are at least 
(21 

second order in and the leading term is given 

by Eq. (10). One realizes that the hadronic decay rate 
should scale as 

h| \KL) X {KL\ i/go ■ (A4) 

For instance and without lost of generality, let us write 
the amplitude of the QCD multipole expansion approach 
to the spin-nonflip two pion hadronic transitions between 
heavy quarkonia: 


Meiei 


. 9 % ^ {<^>F\xk\KL) {KL\xi\<^>i) 
X {h\E^,Er\0), 


(A5) 


where one is assuming factorization of the heavy-quark 
interaction and the production of light hadrons. From 
Eq. (A5), the hadronic decay rate is proportional to 


(fepap’) X (kjai) = [kp {^F\xk\KL)] x [kj {KL\xi\^i )], 

(A6) 


X = (xi + X 2 )l^ is the center of mass position of Q and Q, and 
X denotes xi or X2. 


Initial Final L 

K 

kjr df 

( GeV ) 

= kpo-F 

( fm ) 

ki aj 

( GeV ) 

li 


T ( 4 S ) T ( 1 S ) 1 

1 

1.24 

0.20 

1.27 

0.21 

0.016 

0.017 

1 

2 

1.43 

0.10 

0.74 

0.43 

0.41 

0.89 

1 

3 

1.58 

0.059 

0.47 

0.61 

0.84 

2.62 

1 

4 

1.71 

0.038 

0.33 

0.77 

0.054 

0.21 

1 

5 

1.82 

0.027 

0.25 

0.91 

0.094 

0.43 

1 

6 

1.91 

0.020 

0.19 

1.03 

0.10 

0.54 

1 

7 

1.99 

0.015 

0.15 

1.13 

0.087 

0.50 

1 

8 

2.06 

0.012 

0.13 

1.22 

0.070 

0.43 

1 

9 

2.11 

0.0098 

0.10 

1.30 

0.056 

0.37 

1 

10 

2.16 

0.0081 

0.088 

1.36 

0.045 

0.31 

1 

11 

2.20 

0.0067 

0.075 

1.42 

0.037 

0.26 

T ( 3 S ) T ( 1 S ) 1 

1 

1.24 

0.20 

1.27 

0.44 

0.21 

0.46 

1 

2 

1.43 

0.10 

0.74 

0.66 

0.67 

2.26 

1 

3 

1.58 

0.059 

0.47 

0.85 

0.066 

0.29 

1 

4 

1.71 

0.038 

0.33 

1.02 

0.13 

0.66 

1 

5 

1.82 

0.027 

0.25 

1.16 

0.10 

0.61 

1 

6 

1.91 

0.020 

0.19 

1.28 

0.078 

0.51 

1 

7 

1.99 

0.015 

0.15 

1.38 

0.059 

0.41 

1 

8 

2.06 

0.012 

0.13 

1.48 

0.045 

0.34 

1 

9 

2.11 

0.0098 

0.10 

1.55 

0.035 

0.28 

1 

10 

2.16 

0.0081 

0.088 

1.62 

0.028 

0.23 

1 

11 

2.20 

0.0067 

0.075 

1.68 

0.023 

0.19 

T ( 25 ) T { 15 ) 1 

1 

1.24 

0.20 

1.27 

0.79 

0.44 

1.76 

1 

2 

1.43 

0.10 

0.74 

1.02 

0.16 

0.84 

1 

3 

1.58 

0.059 

0.47 

1.22 

0.13 

0.82 

1 

4 

1.71 

0.038 

0.33 

1.39 

0.089 

0.63 

1 

5 

1.82 

0.027 

0.25 

1.53 

0.061 

0.47 

1 

6 

1.91 

0.020 

0.19 

1.66 

0.044 

0.37 

1 

7 

1.99 

0.015 

0.15 

1.77 

0.032 

0.29 

1 

8 

2.06 

0.012 

0.13 

1.86 

0.025 

0.23 

1 

9 

2.11 

0.0098 

0.10 

1.94 

0.020 

0.19 

1 

10 

2.16 

0.0081 

0.088 

2.01 

0.016 

0.16 

1 

11 

2.20 

0.0067 

0.075 

2.07 

0.013 

0.14 


Table IV. The relative size of the overlap terms that are 
involved in the hadronic transitions T{nS) —>■ T(lS)7r’'"7r“ 
with n = 2, 3, 4. 


and not, as usually considered, to (ka)^ with a the 
average size of the heavy quarkonium states involved in 
the hadronic transition and k the available momentum 
of the process. 

Table IV shows the relative size of the overlap terms 
that appear in Eq. (A6) and are involved in the T(4S') —>• 
T(15)7r+7r“ transition. We compare them with those 
obtained for the T(2S') and T(35') states in which usually 
QCDME approach is considered to work well. One can 
conclude that QCDME approach works at the same level 
for the T(4S') state than for the T(2S) and T(3S) states. 
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